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I In spherically symmetric charged dust, just like in neutral dust, two kinds of 

' singularity may be present: the Big Bang/Crunch (BB/BC) singularity, and shell 

crossings. Quite unlike in neutral dust, the BB/BC singularity may be avoided. 

^ ■ When the charge density pe and the mass-energy density e obey \pe\ < De'= VGe/c^, 

, the conditions that allow the model to avoid the BB/BC singularity necessarily 

I lead to shell crossings. However, when \pe\ — > Dg at the center of symmetry while 
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\pe\ < De elsewhere, both kinds of singularity may be avoided for a sufficiently 
Q^. long period that a body of charged dust may go through the tunnel between the 

^ \ singularities in the maximally extended Reissner - Nordstrom spacetime, to emerge 

into another asymptotically ffat region. An explicit example of such a configuration 
is presented and discussed. It does not contradict any astrophysical constraints. 

X 

: 1 Spherically symmetric dust in electromagnetic field 

For the spherically symmetric spacetimes in the comoving coordinates, the metric can be 
put in the form 

ds^ = e^(*'")dt2 - e^(*'")dr2 - R\t, r) [d^^ + ^\^\^)^^^] . (l.l) 

The case R,r = requires separate consideration. In this case, a solution of the Einstein 
equations for electrically neutral dust was found by Datt [1] and then derived and inter- 
preted by Ruban |2l[3]; the corresponding Einstein-Maxwell equations were systematically 
derived and solved by Ruban |11 O [6] . 

We will discuss the solutions of the Einstein-Maxwell equations in the generic case 
R,r 7^ 0. Most of the material of Sections 2 to 7 has been known for a long time, but a 
few points were overlooked by some of the earlier authors. That material is reviewed here 
in order to introduce the framework and notation, and to fill in the few little gaps. New 
results are reported in Sections 8 and following. 
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29. 



1 



The following problem is considered. In a spherically symmetric charged dust, just 
like in neutral dust, two kinds of singularity may be present: the Big Bang/Crunch 
(BB/BC), and shell crossings. Quite unlike in neutral dust, the BB/BC singularity may 
be avoided. However, when the charge density pe and the mass-energy density e obey 
\pe\ < Df. =^ \/Ge/c^, the conditions that allow the model to avoid the BB/BC singularity 
necessarily lead to shell crossings. (The repulsive term that prevents the BB/BC singu- 
larity in this case is a purely relativistic effect.) There exist initial conditions that allow 
us to avoid both singularities - when either |pe| > everywhere, or |pe| at the 

center of symmetry while |pe| < -De elsewhere. The first case was considered by Ori [7], 
but remained unpublished. In the second case, both kinds of singularity may be avoided 
for a sufficiently long period that a body of charged dust goes through the tunnel between 
the singularities in the maximally extended Reissner - Nordstrom spacetime, to emerge 
into another asymptotically flat region. An explicit example of such a configuration is 
presented and discussed. It has < 0, i.e. it belongs to the recollapsing class that could 
possibly be pulsating for ever. However, shell crossings prevent it from existing for longer 
than one full cycle of collapse and expansion. As shown in the last section, the parameters 
of the object do not violate any conventional astrophysical limitations. 

It is generally expected that shell crossings are just an artefact of the assumption 
that the spatial gradient of pressure is zero. In a more general solution, they would be 
replaced by acoustic waves of a high but finite density that become singularities in the 
limit of spatially homogeneous pressure (in particular, of zero pressure). However, no 
exact solutions of this kind are available, and the singularities in such configurations have 
so far been considered only via existence theorems applied to Einstein's equations with 
various sources, such as the Einstein - Maxwell - Vlasov equations treated in Ref. [8]. 
Typically, the existence theorems predict the properties of solutions over a rather limited 
time-interval. Lacking any more general exact solution, the charged dust models are the 
best device that is available for investigating long-term evolution of charged bodies. 



2 Spherically symmetric dust models with i?,^ ^ 

With spherical symmetry, in the coordinates of (11. ip . the electromagnetic tensor can have 
at most two nonzero components, Ftr (the electric field) and F^<^ (the magnetic field). 
The F^^ 7^ is due to a distribution of magnetic monopoles. Since they are known 
(from experiment) not to exist, and since their inclusion does not lead to any significant 
generalisatioiJl], we will assume the magnetic monopole charge to be zero everywhere. 
The Einstein-Maxwell equations are in this case 
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+ '^^[FJ^F,p^\g^^F,^F^^\^, (2.1) 

F^^, = (47^/c)pen^ (2.2) 
Fv^v,p\ = 0, (2.3) 



^With magnetic monopoles taken into account, the electric charge Q has to be replaced in all formulae 
by \/ Qe^ + Qrn' , whcrc Qe is the electric charge and Qm is the magnetic charge. 
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where e is the mass density of dust, is the velocity field of dust, and pe is the density 
of the electric charge. 
From (12. 2p we find 

= g(r)e-(^+^)/V/^^ (2.4) 

Q„ = (47r/c)pee^/'i^^ (2.5) 

where Q{r) is an arbitrary function; it is the electric charge within the r-surface as (12. 5 p 
shows. With the assumed F23 = 0, eq. (12. 3p is fulfilled identically. 

Now the coordinate components of the Einstein equations with the metric (II. ID and 
the electromagnetic tensor (12. 4p - (12. 5p become 

Goo = ^-" + §|e--Ae-, (2.6) 

Goi = 0, (2.7) 

Gn = -§|J«^ + Ae^ (2.8) 

G22 = G,,/siTi\^) = %^+KR\ (2.9) 



For Gai3 we find 



C4 



_^^c-A _ ^ J ^ ^2.10) 

„A p 2 n u 

A-c ( R^t^ r,R^tt C,tR,t 



V ^ R R J 

G22 = ^ {—4:RR,tt+2RC,t R,t —2RA,t R,t 
-R^Aj-2R'A,u +R^C,tA,t) 
+ — e ^ (ARR^^j. -\-2RC,r Rir — 2RAjj. i?,^ 

+i?^C,j.^ + 2R'^C,rr — R^C,r A, J. ^ . (2.13) 
Equations (12.11) - (12. 3p imply the following: 

{eu^);f, = 0, (2.14) 

eu'';pu^ = —peF'^pu". (2.15) 

Equations (12.141) (the conservation of mass) and (I2.15P (the Lorentz force acting on charges 
in motion and pushing them off geodesic trajectories) are quite general and independent 
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of any symmetry properties of spacetime. Applying fl2.15p to our metric (11.11) we get 



eC„=— (2.16) 
' 27r V ; 



Applying (12.141) to our metric (II. ip we then get 

fei?V/^ = ^Ar_ (2.17) 

where N,r is an arbitrary function of integration. We see that so defined is the energy 
equivalent to the sum of rest masses within the r-surface. From (12. 5p and (12.170 we now 
see that the ratio Q,r /N,r = pe/ (ce) is time-independent. 

Using now (I2.16P and (12. lip we obtain from the equation Gqi = 

2e-^/^i?„. -e-^/M„ R„. = \ ^ . (2.18) 

(It is here that the case R,r = has to be set aside for separate investigation. With R,r = 0, 
the equation Gqi = reduces to R,t C,r = and cannot be used to determine R,r as we 
do below.) Since the expression in parentheses is independent of t, this can be integrated 
with the result 

e-^/^i?„=r(r)-^, (2.19) 
ceR 

where r(r) is an arbitrary function of integration. Using (12.50 and (I2.17P we now find 

PeQ _ QQr 



- QQ,N. (2.20) 

ce A,r 

With this, (12191) becomes 



QQ 

e -'-n,r= i [T) - 

and now (12.160 becomes 



^/'R,r=T(r)-^, (2.21) 



qA/2 

C,r=2^QQ,N. (2.22) 

Using (12.210 and (12.220 to eliminate R,r and C,r from (12.120 and integrating the Gu 
equation we get 

e-On,,^.T^-^,^Jm + m^^^^^hn^ (2.23) 

R R o 

where M(r) is an arbitrary function. Comparing this with the Newtonian equation of 
motion we see that (F^ — l)/2 plays here the role of the energy function E{r). 

The function M (r) is the effective mass that, together with Q, drives the evolution. 
However, as we will see, M(r) is a combination of mass and charge that need not be positive. 
In order to see this, we compare (I2.23P with the Newtonian limit, assuming A = 0. Let 
[x] denote the physical dimension (unit) of x. The dimensions of the quantities appearing 
in (I2:23|l are: [R] = [length], [e^/^dt] = [ds] = [c] x [time], [M] = [G] x [mass]/[c2], 
[A^] = [mass] x [c^], [Q] = [charge] = [^mass] x [length^' ^]/ [time]. The function (F^ — 1) 
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is dimensionless, but, for consistency, must be assumed to have the form 2£/c^, where 
[£] = [velocity^]. In order to find the Newtonian hmit of (12.231) . we muhiply it by and 
let c — > oo. Denoting the Newtonian time by r we obtain 

= 28 + 2Gm{r)/R. (2.24) 

But the Newtonian equation of motion of spherically symmetric charged dust is 



rj = 2S(r) 



GM{t)-P^Q{t) 



/r, (2.25) 



where A4{r) is the total mass within the sphere of radius r, Q{r) is the total charge within 
the same sphere, pe is the charge density and is the mass density. Thus, the M(r) in 
(I2.23P corresponds to the Newtonian [GA4{r) — {pe{r) / Pi^i{r)) Q{r)], and, as announced, 
does not have to be positive. 

To verify the 6*22 equation, we have to find C,t from the Gu equation and A,t from the 
Goi = equation. Substituting these, then finding A,^ from fl2.2ip . using the r-derivative 
of (12.231) to eliminate R,tR,tr, and again using (I2.2ip to eliminate e""^, we obtain 



QQ. 



N 



-^riv„ + (M + gg,^r) 



c4 



0. (2.26) 



One solution of this is Q,n = 0, i.e. a constant total charge. We will mention this simpler 
case later. When Q,j^ 7^ 0, 



c 

The quantity 



^riV„= (M + QQ,^r)„. (2.27) 



M = M + QQ,N r, (2.28) 



an exact analogue of the Newtonian A4{r) of (I2.25p . will appear again in Sec. [31 There we 
will find that Ai is the active gravitational mass. Thus, via (12.271) . F determines by how 
much Ai increases when a unit of rest mass is added to the source, i.e. F is a measure 
of the gravitational mass defect/excess. Solutions with F = are known, this is the Datt 
- Ruban class pLj - |4j. Negative F can also occur (see subsection 18.20.2 of Ref. [B]). 
However, the case F > corresponds to the most ordinary configuration. 

The final equation to take into account is (12.60 . and it reproduces (12.171) . Using (I2.2ip 
to eliminate e^^'^ from (I2.17p . we get e in an equivalent form: 



C^R^R,r V R 

Finally, the Einstein-Maxwell equations for charged dust reduced to a set that defines 
the functions A{t,r), G{t,r) and R{t,r) implicitly. The solution is constructed as follows: 

1. Choose the functions N{r), Q{r) and F(r), and then solve (12.270 to find M(r). In 
fact, since the coordinate r is determined up to the transformations r' = f{r), where / is 
an arbitrary function, we can choose A^, Q or F as the radial coordinate. In the example 
in Sec. [TTl our radial coordinate will be N{r). 
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2. Given these, express e^/^ through R via fl2.2ip . 

3. The set of equations (12:221) - (12:231) then defines e^ and R. In solving (12:22|) 
numerically for C{t,r), an initial condition has to be assumed. Note that so far C is 
not defined uniquely - the coordinate t can be transformed hj t = git'), and then C 
transforms by C = C" — 2 In ((;,(/). Using this, we can require that at the center of 
symmetry C{t,rc) = 0, i.e. for the particle that remains all the time at the center the 
proper time s = t - the time coordinate in spacetime. 

If Q,iy = 0, then Q = const, i.e. Pe = from (12.51) . This case is the neutral dust moving 
in the exterior electric field of a charge concentrated at i? = 0. 

With vanishing charge, eq. (I2.15P reduces to the equation of a geodesic, while eqs. 
(I2.2ip . (I2.22P and (12.230 reduce to those defining the the Lemaitre -Tolman model O I6]. 

3 Matching the charged dust metric to the Reissner— 
Nordstrom metric 

We first transform the R-N metric with A to the appropriate coordinates. For the standard 
curvature coordinates we write {t,Rrn), and we introduce the symbol 

, def 2m 1 . n , , 

h^^l- + + -AR^^\ (3.1) 

^RN rlRN -J 

We demand that the new coordinates (t, r) are still orthogonal, so that gtr = 0, i.e. 

hr^t T,r —-rRRN,tRRN,r = 0. (3.2) 

h 



The function CjiN(t,r) is 



h 

We now solve (13. 2p - (13. 3p for r,^ and calculate 



e^«^ = hr/ - Irrn/. (3.3) 



/ \ — u 2 -'- D 2 G'^RRN,r'^ ^\ 



The component {goo)RN is not fully determined at this point, we will determine it later. 
Since r is defined by the partial differential equation ( 13. 2p . it still involves an arbitrary 
function of one variable. 

We wish to match the charged dust metric of Sec. 1 to the R-N solution given above 
across a hypersurface r = r^,. This requires that the induced 3-metric and the second 
fundamental form of this hypersurface are the same for both spacetime metrics. Continuity 
of the 3-metric requires that 

^c{t,r,) ^ ^CaM{t,n)^ ^^^^^^^ ^ RRN{t,n). (3.5) 

The transformations that keep the metric diagonal are still allowed. Transforming t by 
t' = J e"(*Mt, where a = CR]\f{t,rh) — C(t,rh) we fulfil the first of (13. 5p . while g[^ and gQ^ 
are not changed. 
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On the surface r = r^, RjiN(t,rb) must be the same function of t as R{t,rb). Conse- 
quently, R,t (t,ri,) must be the same in both metrics, and so Rji]\f(t,ri,) must obey 



^RRN,t{t,n) = T'{n)-l + 



RRN{t,rb) C^RRN^{t,n) 



■^ARn^^t^n). (3.6) 



The unit vector normal to the hypersurface r = Vb has components 

= (0, e-^/^ 0, 0) ^ (o, (r - ^) / R,, , 0, 
for the interior metric, and, from (13 .4^ . 



(3.7) 



.yg«=(° /''^;"""^ Mo,0| (3.8) 

^RN,r 

for the R-N solution. From the continuity of the second fundamental form we have 
R,rXX^^^ = RRN,rX^RN\r=r, ^nd (e^),.^^^^^^ = (e^«- ) ^Ijiv | • The first con- 
dition says 

(/. + e-^-i?«^/)_^ = (r - , (3.9) 

which ensures the continuity of gn = —e^ across r = r^, even though we have not required 
this. Substituting for Rrn^i^ from (13.61) and for h from (13.11) . then comparing the coefficients 
we obtain 

e = ^Q(r,), m = (M + QQ,^r)^^^^. (3.10) 

The continuity of the second fundamental form imposes one last condition, on C,r. 
Using (1223), fl^A and ([S3D, the condition is 



T-QQ,n/R 

^RN,r 



R 



RN,r 



= 2^. (3.11) 

r=rt, 



We have no expression yet for Crn^t-, and we will find it from the field equations now. We 
know that e^ is continuous at r = r^, so we can use (I2.2ip for AR]s[{t,rb). Substitute this 
in (12. lip , and take the equation Gqi = at r = r;,. The result is 



^ /, N '^RRN,rQQ,N 



RHT-QQ,M /R) 



(3.12) 

r=r5 



and it shows that ( 13. lip is fulfilled. 

Thus, (I3.10p are the only limitations imposed on the charged dust metric by the match- 
ing conditions. This matching was first discussed by Vickers [9]. The second of (I3.10p 
reveals the connection between the active gravitational mass (in this case m) and the 
effective mass (compare (I2.28P ). 
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4 Prevention of the Big Crunch singularity by electric 
charge 

In the present section we will deal only with the Big Bang/Big Crunch singularities, as- 
suming A = 0. The shell crossings will be discussed separately in Sec. [71 In the following 
we shall denote E{r) = {T^ - 1) /2. 

The presence or absence of a singularity is detected by investigating the roots of the 
right-hand side of (12.231) . which, for this purpose, is more conveniently written as 

e-^R^R,t^ = 2E{r)R^ + 2M{r)R + {Q,n^ - G/c^) = W{R). (4.1) 

At each root of W{R), the sign of R,t changes, and evolution is possible only in those 
regions where W{R) > 0. The following cases occur 

(a) When E <0, W{R) has roots only if 

With no roots, W{R) would be negative at all R, so (14.21) is the condition for the existence 
of a solution of (14. ip . The roots are 

R± = -^± ^^M^ - 2EQ^ {Q,^' - G/c'), (4.3) 

and W{R) > between them. Nonsingular solutions will exist when both R± > (with 
R± < 0, no solution exists at all, and with R^R^ < 0, i? = is in the allowed range.) 
This is equivalent to 

Q,N^ < Gjc^ and M > 0. (4.4) 

We will interpret this condition later on in this section. 

If there is equality in (g^D, then W{R) <^ for all i? ^ i?_ = and iy(i?±) = 0. 
Then R = R± and the model is static. If, in addition, Q,]y (r^) = (meaning Peiji) = 0) 
and r(rft) = 0, then E = —1/2, and in this case the exterior R-N metric is the extreme 
one, with = m^. 

With (14. 4p fulfilled, R oscillates between a nonzero minimum and a maximum. 

(b) When E = 0, singularity is avoided if and only if M > and Q,n'^ < G/c^. Collapse 
is then halted and reversed once and for all. 

(c) When E > W{R) > either everywhere (if there are no roots) or beyond the 
roots. There will be no roots when < 2EQ^ {Q,n'^— G/c'^), in which case W{R) > 
for all R including R = 0, and the model can run into the singularity. Thus (14. 2 p is here one 
of the necessary conditions for the existence of nonsingular solutions. With (14.20 fulfilled, 
W{R) has two roots, and at least one of them has to be positive if singularity is to be 
avoided. With M > 0, we have R- < always and R+ > if and only if Q,n'^ < G/c^. 
With M < 0, i?+ > and R- < i?+ always, so nonsingular solutions exist with no further 
conditions, provided R > i?+ initially. Collapse is then halted and reversed as in case (b). 
The bounce with M < is nonrelativistic, since it occurs also in Newton's theory, under 
the same conditions. 

Now we will interpret the condition (14. 4p . The inequality Q,n'^ < G/c^ translates into 
\pe\ < VGe/c, which means that the absolute value of the charge density is sufficiently small 
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(but nonzero) compared to the mass density. This kind of bounce is purely relativistic, and 
it does not occur in the Newtonian hmit: with M > 0, i? = is always in the allowed range 
of Newtonian solutions. The interpretation of the relativistic bounce in the Newtonian 
terms is this: as seen from fl2.23p . the charges provide a correction to the effective mass 
M, so that it becomes M = M + (1/2)Q2 (g,jv^ - G/c^) / R. This correction is negative at 
small charge density (when Q^n"^ < G/c^), so it weakens gravitation, thus helping the dust 
to bounce. However, at large charge density {Q,n'^ > G/c^), charges enhance the effective 
mass and thus oppose bounce. (A similar phenomenon is encountered in the motion of 
particles in the Reissner-Nordstrom spacetime, where an electric charge in the source of 
the gravitational field creates effective antigravitation, provided the charge is small enough 
compared to mass.) Nevertheless, in this last case, the Newtonian electrostatic repulsion 
can prevail, provided M < at the same time. 

With Q,N= the BB/BC singularity is avoided in every case when a solution exists. 
Thus, for neutral dust moving in an exterior electric field, the BB/BC singularity never 
occurs. This was first found by Shikin [TO]. This is a purely relativistic effect. 

The above implies that with (14.41) fulfilled a solution of (12.231) . for which R ^ initially, 
does not go down to 0. However, if the charged dust occupies a volume around the center of 
symmetry R= 0, then, at any time, there are dust particles with all values of R, including 
R = 0. (We will find in the next section the conditions to be obeyed in order that the 
center is nonsingular.) Thus, the inner turning points given by (14.31) will exist arbitrarily 
close to the center. We will use this remark in Sec. [3 

If EM < and Q,n'^ = G/c^, then (H?T|) has the time-independent solution R = -M/E. 
In this case, the electrostatic repulsion just balances the gravitational attraction and the 
whole configuration is static - but unstable. When M > > E, the perturbation can only 
be toward smaller R, and it will send the dust into collapse that will terminate at i? = 0. 
When M < < E, the perturbation can only be toward larger R, and it will send the dust 
into infinite expansion. Another time-independent solution is E = = M, Q,n'^ = G/c^; 
in this case GT /c^ = D = const and = DAi . 

The surface of the charged sphere obeys the equation of radial motion of a charged 
particle in the Reissner-Nordstrom spacetime. For such a particle, if the ratio of its charge 
q to its mass obeys {q/ fJ^y < 1, then the reversal of fall to escape can occur only inside the 
inner R-N horizon, at i? < r_ = m — \/m^ — e^. Thus, the surface of a collapsing sphere 
must continue to collapse until it crosses the inner horizon, and can bounce at -R < r_. 
Then, however, it cannot re-expand back into the same spacetime region from which it 
collapsed, as this would require motion backward in time, as seen from the compactified 
spacetime diagram of the maximally extended R-N solution [6j. The surface would thus 
continue through the tunnel between the singularities and re-expand into another copy 
of the asymptotically flat region. This possibility is interesting, since it shows that the 
maximally extended R-N spacetime, usually interpreted as an abstract geometric curiosity, 
may in this case become an astrophysical reality. 

The bounce at small charge density {Q,n'^ < G/c^) is more interesting physically, since 
the real Universe has no detectable net charge, so only small localized charges could exist 
in it. We saw that an arbitrarily small uncompensated charge can prevent the BB/BC 
singularity. Unfortunately, Ori [HI [12] proved that if Q^n"^ < G/c^ holds throughout the 
volume, then a shell crossing is unavoidable, and it will block the passage through the 
tunnel. We will derive this result in Sec. [71 Thus, a nonsingular bounce through the RN 
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tunnel is possible only if Q^n'^ > G /c^ everywhere or if Q^n^ G/c^ at the center, while 
Q^n'^ < G/c^ elsewhere. 

5 Charged dust in curvature and mass-curvature co- 
ordinates 



It is instructive to transform the metric given by (11. II) with (12.211) - (12.231) to such co- 
ordinates in which the function R[t, r) is the radial coordinate. We note that R,r dr = 
d-R — R,t dt, and we take t to be a function of the new coordinates: t = f{T, R). Thus 



R,r dr 



dR-R,t{f,rdT + f,RdR), 

dt = f,,dT + f,RdR, 



and the new metric components, using (12.231) . are found to be 



9tt 
dRR 



e^f,,f,nA + f,,R,t 



2 ' 



{T-QQ,^/R) 
e^f,R^A-l + 2f,nR, 



{T-QQ^n/RY 

+ 



R c^R^ 3 



(5.1) 

(5.2) 
(5.3) 

(5.4) 

(5.5) 



where A4 is defined in (I2.28p . Note, from ( 15. ip . that the transformation is different for 
collapsing dust {R,t< 0) and for expanding dust {R,t> 0). The transformation from 
(r, i?) to (t,r), inverse to (15.11) . is analogous to introducing, in the Reissner-Nordstom 
region r > r+ = m + y/m'^ — e^, coordinates comoving with the congruence of charged 
particles that are radially collapsing or expanding, respectively. The extension is to the 
future or to the past, respectively. 

This can now be specialised in two ways. One possibility is to choose the proper 
curvature coordinates, in which QrR = 0. This representation of the Vickers metric has not 
found any application in the literature so far, but it is instructive - see Ref . ^ . The other 
possible specialisation of (15. 2p - (15.41) is to choose the defined in (I2.28P as the new r 
coordinate. These mass-curvature coordinates were first introduced by Ori [11]. The 
surfaces A4 = const are timelike, so none of the coordinates is time and the metric in these 
coordinates cannot be diagonal. Since A4, Q, T and depend only on r in the original 
coordinates, we have Q = Q{M.), T = r(A1) and N = N{M.). The Jacobi matrices of the 
transformations {t,r) ^ {A4,R) are 



d{t,r) 
d{M,R) 



f\M , f\R 

r,M , r,R 



d{M,R) 
d{t,r) 



M,r 
R)t 1 Rir 



(5.6) 



These matrices must be inverse to each other, hence 



R, 



M ■ 



R,tM,r' 



f,R=l/R,t, r,M=l/M, 



r,R=0. 



(5.7) 
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In the coordinates {x^,x^) '= {A4,R) the velocity field still has only one contravariant 
component: 



„H.±,/r-i + ^ + ^^(^^f^-iA«'^ (5.8) 

(+ for expansion, — for collapse). We define the auxiliary quantities 

u = T- QQ,N IR, (5.9) 

{f'^lu)f.M = F{M,R), (5.10) 
and, using ([ESD, ([EID and ([Z^SD, we get in ([ED - CT : 

gMM = F^^, gMR = Fu/u'', gRR = l/{u'')\ (5.11) 

while eqs. flETl) simplify to dr = dM/M,r and dt = {-R,rdM/M,r+dR) / R,f The 
function F{A4, R) is to be found from the field equations. 

Using (15. 7p . (12. 4p and (IS.lOp we find that the only nonvanishing components of the 
electromagnetic tensor in the {Ai,R) coordinates are 

r^MR _ TpRM _ Q p _ p _ ^K^0\ 

^ - ~^ - ~ -J^RM - (5.12; 

Further, using ([MD, O, ([MID, (EiSTD, (12:231) and ([EH]) we find for the charge 

density and energy-density 

47rpe _ Q.M _ 2 

Recall that the (A^, -R)-coordinates cover only such a region where R,t has a constant 
sign. As seen from (15. 7p and (I5.10p . F changes sign where does, and so does = 
e-^/^R,t. Thus, Fu^ preserves its sign when collapse turns to expansion and vice versa. 
This observation will be useful when calculating F later. 

Now F{Ai, R) is the only unknown function. We obtain further 

um = uF, ur = (5.14) 

{Fu^f Fu^ ^ ' 

In the (tVI , R) coordinates eq. (12.270 reads 

(GT/c^)iV,^ = 1. (5.16) 
The function F can be found from (I2.15p . Using (I2.20p . it is found to be: 



Using (15. 9p and (15.161) we transform (I5.17P to 

1 r„ 1 



(5.17) 



,R'- 



r,A4 +- 
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{Q,N + QQiNN ) 



(5.1^ 



This coincides, except for notation, with Ori's (1990) result. As Ori [llj stressed, eqs. 
(15. lip . (15. 8p and (I5.18P determine the metric exphcitly, in contrast to the representation 
by Vickers used in Sec. [21 where the Einstein-Maxwell equations were reduced to a set of 
two differential equations. However, we lost the information about the time-dependence 
of R. Points of the spacetime are now identified by the values of Ai and i? - by specifying 
the pair {Ai,R) we say: this is the point in which the shell containing the mass Ai has 
the radius R. However, we have no means of saying, without recourse to the comoving 
coordinates, how much coordinate time it has taken the shell to expand from the minimal 
size R = -Rmin to the current R. In consequence of this, the {Ai,R) coordinates do not 
allow us to see whether the minimal size (which is a singularity or a nonsingular bounce) 
was achieved by all shells simultaneously with respect to the time coordinate t or not. This 
information is crucial for considering shell crossings, as we will see in Sec. [71 

With A = 0, the integral of (15.180 is elementary, but requires separate treatment of 
various subcases. The full list of results is given in Ori's paper. 

Now the Einstein-Maxwell equations are all fulfilled. 

Note from (I5.1ip and (I5.13P that the metric and the mass density are insensitive to the 
sign of u^. However, as explained in the remark after (15. 5p . > and < correspond 
to different maps with different domains. Thus, integrating (I5.18P from RitoR> Ri with 

> 0, we integrate forward in time, while calculating the same integral with m-^ < we 
integrate backward in time. 



6 Regularity conditions at the center 

Just as in the L-T model, the set i? = in charged dust consists of the Big Bang/ Crunch 
singularity (which we showed to be avoidable) and of the center of symmetry, which may 
or may not be singular. We will now derive the conditions for the absence of the central 
singularity. We assume no magnetic charges. 

Let r = Tc correspond to the center of symmetry. From (II. ip . (12.210 and (12.290 we see 
that N{r) = e^/^d3X = Air /J' e{r')R^{t, r')dr', where V is a sphere centered at r = Tc in 
a t = const space. Thus, if e has no singularity of the type of the Dirac delta at the center, 
N{r) must obey N{rc) = 0. Similarly, eqs. (11.10 and (12.50 show that if there is no delta-type 
singularity of pe at the center, then the electric charge must obey Q{rc) = 0. With both e 
and Pe being nonsingular at Tc, and e{t,rc) > 0, the ratio Pe{rc) / e{rc) is nonsingular, and 
(12.50 with (I2.17P show that lim^^r-c Q/^ = ^i^r^r^ Q^n = Pe{rc)/[ce{rc)] is finite (possibly 
zero). Then, (15.160 implies that r(rc) 7^ and lim^^r-c N/ M. = limr^.^^ N,_m = GF/c^, and 
this, together with (rc) = 0, implies Ai (r^) = 0. Then, from (I2.28p . also M (r^) = 0. 

Since i?(t,r,) = and A^(r,) = 0, we find from (K29\\ : 

r r ^R^R-,T y 3 / QQ-,n\ 

hm — = hm = hm 1 . (o-l) 

r^Tc N r^rc N,r r-*rc Ane \ R J 

This limit will be finite if \im.r-,rXQQiN / R) < 00. We already know that Q,n (fc) == 5o = 
const and thus \imj^_,Q Q / N = %. Thus, \\mr^r^Q / Ai = go and limr^rAQQ^N I R) < 00 
if \im.r^rc Rj AV = const, where 7 < 1. Then, (16.10 imposes the further condition that, in 
the neighbourhood of Vc 

lim R/M^^^ = Pit). (6.2) 



12 



We assume P(t) ^ 0, since /? = implies, via (16. ip . one of two unphysical situations: (I) 
lim^^j.^ e = oo, i.e. a permanent central singularity, or (II) lim^^rc (r — QQ^n /-R) = 0, 
which leads to e(rc) being independent of time - a pathological situation in an expanding 
or contracting configuration. 

Since i?(t, r^) = at all times, we have R^t {t, f^c) = 0, and, from the above, limj-^^c 
R,t I M^l"^ = P,t- AH other terms in (I2.23P except (F^ — 1) vanish at r = Vc, and so we 
must have 

lim r^(r) = 1 ^ lim E{r) = 0. (6.3) 

r—>rc r—>rc 

From the limiting behaviour of the functions in (12.230 at r ^ r^, we conclude that 

lim lEjM^I^ = lim (T'^{r) - l) /A^^/^ = const, (6.4) 

r— >rc r— +r'c 

and this constant may be zero. Note that with F < this means F(rc) = —1. A central 
point where F < is the "second center of symmetry", in those models that have it. It 
is the antipodal point, in the spherical space, to the ordinary center. Having reached the 
antipodal point, we have added as much mass to the space as it can contain, and adding 
new mass is not possible. The condition F(rc) = —1 must then be understood as follows: 
the active gravitational mass increases when we take away the rest mass from the object. 



7 Shell crossings in charged dust 

As is seen from (IS.lOp and (15.70 . F = is equivalent to R,r = 0, so F = is a locus of shell 
crossing. Then, from (15.130 we see that TFu^ must be negative for the density of dust to 
be positive. Since = dR/ds < during collapse, TF must then be positive. 

We will investigate the occurrence of shell crossings in a configuration that avoids the 
Big Bang/Crunch singularities. Let us write the solution of (I5.18P as follows 



F = I{M,R,Ri)+g{M), (7.1) 



where 



R 



1 r„ 1 



-^{Q,n +QQ,nn) 



dx, (7.2) 



Ri is the initial value of R and g{Ai) is an arbitrary function - the value of F{Ai,R) 
aX R = Ri. In following the collapse of the central region forward in time, we eventually 
come close to the bounce point, where R < Ri. The integral becomes unbounded as R 
approaches the bounce value (where ^ 0), and so, with R < Ri, the first term in (17. ip 
is positive when the expression in curly braces is negative. 

We saw in Sec. Hlthat there are dust particles with all values of R, including R = 0. At 
the turning point 0. The integrand is of the form (ax^ + bx + c) '^^^ (02^'^ + OsX^) , 

and the trinomial has real zeros, so the integral is unbounded, which shows that = 
is a coordinate singularity!! From (17. 2p it is seen that as we get near to {R = 0, = 0}, 



^Incidentally, the set where F = is also a coordinate singularity. This set is the "neck", known from 
studies of the Lemaitre-Tolman model, see Ref. [13] . 
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the coefficient [1 / (xT)] will become unbounded. We know from the regularity conditions 
that limj.^rc = 1, lim^^rc Q = 0, lim,.^rc Q,n = const < oo and limr^,,^ Q/R = 0. Thus, 
as long as Q,n'^ < G/c^, the term containing l/(xr) will dominate in the vicinity of i? = 
and will determine the sign of the infinity in F. Now it turns out that the sign of F will 
necessarily change to opposite during collapse: If F > and g{M) > 0, then F > at 
R = Ri, but F — oo as i? — > 0. If F < 0, then F < can be achieved aX R = Ri by 
the choice g{Ai) < 0, but F +oo as i? ^ 0. This means, there is necessarily a shell 
crossing somewhere at i? > if Q,n'^ < G/c^ holds all the way down to = 0. This is 
the theorem proven by Ori pTl WI\ . 

The infinity in F could be avoided if the term in curly braces in (17.21) were zero at the 
same x, at which = 0. The zero of is given by (14.31) . it is i? = In that case 
F is finite at i? = and, by (I5.13p . e becomes infinite, i.e. R = R+ becomes a true 
curvature singularity. Thus, also in this case, the charged dust cannot tunnel through the 
Reissner-Nordstrom throat. 

The only situations in which both BB/BC and shell crossing singularities could possibly 
be avoided are these: 

1. When limr^rcQ^N^ = G/c^ - then, because of ]Sm.r-,r^Q / R = 0, the term F,_v! in 
(17.21) has a chance to outbalance the other one and secure the right sign of F everywhere. 

2. When Q,7v^ > G'/c^ F > and M < 0. Then, as seen from (gSD, both /2± > and 
a nonsingular bounce at with no shell crossings, is possible. 

Ori (private communication) found an example of a fully nonsingular bounce with 
Q^n"^ > G/c^, but has never published this result. An example of a nonsingular bounce 
with Q,n'^ < G/c^ will be given in Sec. [TTl 

8 The evolution of charged dust in mass- curvature 
coordinates. 

In the mass-curvature coordinates, with A = 0, it is easy to solve the evolution equation 
di?/ds = u^. By this opportunity, we can again identify, by another method, all those 
solutions that avoid the BB/BC singularity. 

We take £ = +1 for expansion and £ = —1 for collapse jf] we denote 

=^ F^ - 1 $ = 91S^'^^ ~ 1 N 



2E 



and obtain from (15.81) with A = 



d R 

— = iJ2E + 2M/R + 2E$/i?2. (8.2) 
ds 

The integral is different for each sign of E, and for E > there are separate subcases 
depending on the values of $ and M; $ < corresponds to Q,n'^ < G/c^. With $ = the 
solutions are of the same form as those in the L-T model, although, if Q 7^ 0, the effective 
mass M still contains a contribution from the charges. With 7^ it is most convenient 
to represent the solution R{s) by parametric formulae {R^u), s^u)}. 



•^We refer here to the initial instant of evolution. In the nonsingular models, collapse will reverse into 
expansion during the evolution. 
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Figure 1: Schematic graphs of the solutions (|8.3p . with different values of M and $. Only those parts 
of the curves that lie above the s axis describe physical situations. From bottom to top (counted at 
the left edge), the curves correspond to: 1. {I2 = — 1,M < 0,$ > 0}. This model has a finite time of 
existence. 2. {^2 = +1, M > 0, $ < 0}. 3. {^2 = +1, M > 0, $ = 0}. This is the uncharged (Lemaitre - 
Tolman) model, included for comparison. 4. {I2 — +1,M > 0,$ > 0}. The collapsing branch (left) ends 
in a singularity at i? = 0, the expanding branch (right) starts at the singularity at i? = 0. Models 1, 3 
and 4 are the only ones with singularities. 5. {I2 = +1, M < 0, $ > 0}. 6. {ii = +1, M < 0, $ = 0}. 7. 
{€2 = < 0, $ < 0}. The instant of time-symmetry is in general different on different Af-shells. It 

was made the same in the figure only for better readability. 



With > and $ < MV {\E'^), the solution is 

R = 
s-ssiM) = 



M2 ^ , M 

— <P coshcj , 

4^2 2E' 




- <I> smh uj — uj 

4^2 2E 



.3) 



where u is the parameter, £2 = =tl and sb{-M) is an arbitrary function of integration. 
In the L-T hmit Q = 0, this function becomes the bang-time. In the bang-free models, 
sb{-M.) represents the instant at which R achieves a minimum or a maximum. 

The solution (18. 3p has no BB/BC singularity when £2 = +1 and either M < or $ < 0. 
Schematic graphs of the solutions are shown in Fig. [H 



With E > and $ = M2/ (4^2), the solution is 



s-SsiM)^-^R-^ln 



R + 



M 
2E 



Rq 



14) 



This one will have no BB/BC singularity only if M < 0. Schematic graphs of these 
solutions are shown in Fig. [2l 
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Figure 2: Schematic graphs of the sohitions (|8.4p . with £ = +1 and different vahies of M. Curve I is 
the nonsingular solution with M < 0, curve II is a sohition with M > and a singularity at finite s. In 
solution I R tends to zero asymptotically as s ^ — oo. 



This model always runs into a BB/BC singularity at sinlicij = M/ \J 4i?2$ — M^, inde- 
pendently of the sign of M. It has no uncharged limit $ = 0. Schematic graphs of the 
solutions are shown in Fig. [31 

With £■ = 0, we have to separately consider the subcase M = 0. Unlike in neutral dust, 
this subcase is not vacuum (see fl2.29p ): it corresponds to the case when the electrostatic 
repulsion among the dust particles exactly balances the gravitational attraction, so that 
the effective mass is zero. This subcase exists with all signs of -E, but does not require 
separate treatment when E (Such cases exist also in Newton's theory. The Newtonian 
solutions with M = and > are expanding or collapsing with constant velocity, 
R = ±y/2E {t - to) + Ro, those with M = = E are static.) 

When E = and M ^ 0, the solution is 



s - ssiM) = -—^2MR + g2 (g,^2 _ ^/^4) ^mR - {Q,^' - G/c')] , (8.6) 



and it will not run into a BB/BC singularity only if M > and Q,n'^ < G/c^. Schematic 
graphs of the solutions are shown in Fig. HI 

When E = = M, the solution exists only when Q,n'^ > G/c^, and it is 



With E > and $ > M7 (4^^), the solution is 




(8.5) 




(8.7) 
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Figure 3: Schematic graphs of the solutions (|8.5p . with £ — +1 and different values of M and <i>. Every 
solution has a singularity; only the slope of the curve at the singularity and the instant of the singularity 
depend on M and <&. 




Figure 4: Schematic graphs of the solutions (|8.6p . with £ — ±1 for curve I and £ = +1 for the other 
two curves, and with different values of M and Q,n- Curve I is the nonsingular solution with M > 
and Q,N^ < Gjc^, curve II corresponds to M < and Q,Ar^ > G/c**, curve III corresponds to M > 
and Q,Ar^ > G/c^. Model II achieves maximal expansion at i? = —Q^ {Q,n^ — G/c^) /(2-M") and then 
recoUapses. 
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Figure 5: Schematic graphs of the solutions (|8.8p . with £ = +1 and with different values of M and 
<&. Only those parts of the curves that lie above the s-axis describe physical situations. The values of 
the parameters are as follows: Curve I ~ Af < and $ < 0; Curve II - M > and <i> > 0; Curve III - 
M > and $ = (this is the Lemaitre - Tolman evolution, shown for comparison); Curve IV - M > 
and 3> < 0. Curve II is the nonsingular model, all the other models begin and end with a singularity, and 
the derivative dR/ds is always infinite at the singularity. Different Af-shells have different periods in s. 
The periods were made equal only to improve the readability of the picture. 



This one always has a singularity. A plot has a similar shape as curve III in Fig. HI 
When E < 0, a solution exists only with $ < M^/ (4i?^), and it is 



R 



s - sb{M) 

This will avoid a BB/BC singularity only if M > and $ > 0. Schematic graphs of 
the solutions are shown in Fig. O This is a periodic solution, and the period is T = 
27rM/(-2E)3/2. 

Note that if the functions M(r), E{r) and Q{r) are chosen in accordance with the regu- 
larity conditions of Sec. 6, then the regular behaviour of R at the center, lim^^r.^ R/ Ai^^'^ = 
(3{t), automatically follows in eqs. ([HJD - dUD and dH]), but not in flHTD . 



M 
~2E ~ 

I 




— $ sin ijj 
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9 Transition from collapse to expansion by nonsingu- 
lar bounce 



We will now investigate transferring the solution (17. ip - fl7.2p from the collapsing sector 
to the expanding sector. This is a nontrivial problem because, as we have seen, the 
bounce point along each worldline, where = 0, is a singularity of the mass-curvature 
coordinates, and we deduced flT.ip - fl7.2p having < in mind. On the other hand, all 
equations up to (15.181) apply independently of the sign of u^. 

Let us rewrite eq. (I2.23P treating (A^,_R) as the independent variables and t as the 
unknown function. Proceeding as in eqs. (15. 6p - (15.70 . we find that (I2.23P is equivalent to 

dt e-CiM,R)/2 

(9.1) 



dR 



P - 1 + 2M/R + Q2 (Q,^2 _ ^/^4) _ (i/3)Ai?2 



where + corresponds to expansion, — to collapse, and (F, M, Q) are now functions of Ai. 
The solution of (19. ip is 

t = ±JiM,R,Ro)+tBiM), (9.2) 
where tB{Ai) is an arbitrary function and 

, , fR -C(M,w)/2 i 

J{M,R,Ro)= / = (9.3) 



Ro 



- 1 + 2M/W + g2 (g,^2 _ ^/^4) _ (1/3)A^2 



The constant Rq is an arbitrarily chosen initial point of integration, and tsiAi) is the 
integration "constant" , an analogue of the Lemaitre - Tolman bang-time function. Here, 
tBi-M) = t\R=Ro- Since R is the only quantity that may depend on t in the comoving 
coordinates, we see that it depends on t always through the combination {t — ts)- 

Note that eqs. (19.11) - (19.31) apply only throughout a single expansion phase or a single 
collapse phase. Thus, when we want to consider the transfer from collapse to expansion in 
those models in which it is nonsingular, we have to choose Rq = Rmm, the minimal value 
achieved by R for a given A4, and thus t = Ib will correspond to the bounce instant. 

Consider now a pair of solutions (19. 2p - (19.30 . the first one collapsing, the other one 
expanding, both with the same function tB{M.)i and consider the two solutions at the 
same values of (tM, i?, Rq). In accordance with the previous paragraph, we have to choose 
tB^M.) = t{M., Ruan) for both. Let the collapsing solution go through these values of 
(A1, i?, Rq) at if: = tc, and the expanding one at t = te- We have tc = —J{Ai, R, Rq) + tB 
for the collapsing solution, and te = +J{Ai,R,RQ) + tB for the expanding one. Thus, 
tc + te = '^tB, i.e. tB^J^) is the midpoint of the interval [tc,^e]- In general, it thus depends 
on A^, i.e. is different for each mass-shell. The case t_B = const, i.e. when the midpoint 
has the same value for all masses, corresponds to the evolution that is time-symmetric 
with respect to t = tB- 

Now consider eq. (15.101) (recall: the / in (I5.10p is the comoving time t). Let us consider 
a pair of curves (Cc, Cg) in the {R, ^A) surface given by the same equation R = Rce{M.), 
one curve in the collapsing region, the other one in the expanding region. Write eq. (I5.10p 
first for the collapsing solution, along Cc'- 

qC/2 

^ [-J,MiM,R,RQ)+tB,MiM)] = UM,R,Ri)+gciM), (9.4) 
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where the functions I and g are those defined in (17.11) - (17. 2p . and the subscript c refers 
to collapse. Write the same equation for the expanding solution, along Cg: 



[+J,M {M, R, Ro) + tB,Mi^)] = UM, R, Ri) + geiM). 



(9.5) 



T-QQ,N /R 



The Jg in (19. 5p differs from the Ic in (19. 4p only by the sign of u^, otherwise all quantities 
in /g are the same as their counterparts in Jg. Thus, along the chosen pair (Cg, Cg), we 
have /g + /g = 0. Knowing this, and adding eqs. (19. 4p and (19.50 . we obtain 



This is the equation that relates g{M.) for the collapsing solution in (17. ip - (I7.2p to the 
corresponding function in the expanding solution. For a time-symmetric evolution, and 
only then, we thus get ge = —gc- Note that if the l.h.s. of (19.60 is positive and large, while 
gc is positive and smaller than the l.h.s., then ge cannot be negative. This observation will 
prove important below, in analysing the existence/avoidance of shell crossings. 

We recall (see the comment after (I5.13P ) that F > where Tu^ < and F < where 
Tu^ > 0. As the bounce point is approached from the > sector in the central region, R 
eventually becomes smaller than Ri, F becomes unbounded, and it must tend to — oo. We 
concluded before that the expression in curly braces must be negative if shell crossings are 
to be avoided in collapse. Thus, the first term in (17.10 - (17. 2p gives the correct behaviour 
of F in approaching the bounce point from the > sector. 

The function F is thus discontinuous at the bounce set, it tends to +oo when this set 
is approached from the < sector, and to — oo when it is approached from the > 
sector. But, as seen from eqs. (15. lip , the metric components depend on F^, F/u^ and 
(u^) , and so they do not jump from one infinity to the other. Then, since Fu^ is finite at 
bounce and does not change sign, we see from (15.130 that the mass density is continuous 
across the bounce set. 

Let us note that only the integral term / in (17.10 simply changes sign when we go over 
from collapse to expansion. The function g{Ai) transforms by (19.60 . and it changes to its 
negative only for time-symmetric evolution. Thus, if we have guaranteed that there are no 
shell crossings before the bounce by choosing such g that gives I + g > in the collapse 
phase, but the evolution is not time-symmetric (tB,M 0); then I + g < will in general 
fail to hold in the expansion phase, and shell crossings will appear after the bounce. We 
shall come back to this question later. 

10 Nonsingular bounce of a weakly charged dust 

From now on we assume Q,n'^ < G/c^ (weakly charged dust). 

The cases E > and E < have to be considered separately. In the first case, if 
a totally nonsingular model existed, it would be first collapsing from an arbitrary initial 
density, then it would go through just one bounce, and re-expand so that each non-central 
shell achieves an infinite radius. In the second case, each shell would oscillate between the 
two turning values of R. At each turning value, = and \F\ — > oo, and we have to 



2tB,Mi^)=9c{M)+ge{M). 



(9.6) 



T~QQ,m/R 
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choose the arbitrary functions in (17.11) - (17.21) so that the infinity has always the right sign 
(e.g. F +00 when < and F > before the turning point). Then, in principle we 
can secure the right sign in two ways in each case: by imposing conditions on r,_A/( or on 

Fi ''^ 1 - (cVg) (g,jv' + QQ.NN ) , (10.1) 

depending on which one dominates in a given situation. 

We first observe that a weakly charged nonsingular dust distribution with E > 
cannot exist. We assume that the regularity conditions at the center hold, and that 
Q,n'^ — > G/c*. Throughout the proof, and then in deducing the conditions for no shell 

crossings in an ii^ < model, we assume < (collapse), just for definiteness. In the 
end, it will be easy to see that the same conditions must apply in the expansion phase. 

Suppose, for the beginning, that E = Q and consider collapse. Then F = \/2E + 1 = ±1 
is constant throughout the spacetime, and T^j^ = 0. Thus, Fi must determine the sign of 
F at the bounce point. To secure the right sign also in the central region it must obey, 
in a vicinity of the center, Fi < 0, independently of the sign of F (because TF must have 
the same sign as long as does not change sign). In a vicinity of the center the quantity 
(1 — c^Q,]\i'^/G) is positive by assumption, while Q,n'^ = G/c^ at the center. We have to 
consider two cases: 

1. Assume Q,n= +VG/c'^ at the center. This means that Q,n> in some vicinity 
of the center. Since, by the regularity conditions, Q = at the center, we have Q > 
around the center. Then the condition Fi < means Q,nn> {G/c'^ — Q,n'^) /Q > 0, 
which implies that Q,n in the vicinity of = is larger than at = 0, i.e. larger than 
■\/G/c^. This is a contradiction with the assumption. 

2. Assume Q^n = —VG/c^ at the center. Then Q,n < in some vicinity of the center, 
and from the regularity conditions Q < in the same vicinity. Then Fi < implies 
Q,NN < {G/c^ — Q^n"^) /Q < 0, which means Q,n < —^/G/c^ in a neighbourhood of the 
center, in contradiction to the assumption. 

Thus, E = leads to a contradiction and is thereby excluded. 

Suppose then that E > and that Fi/(_RF) becomes negligible as x — *• 0, so that the 
sign of F,;k determines the sign of infinity of F at bounce. We still consider collapse. Since, 
with F > and < 0, F must be positive, T,m < is the right choice. The regularity 
conditions require that F^ = 1 at the center. If F > 0, then F = +1 at the center, and then 
F,;v4 < means < F < 1 in the vicinity of the center. But 2E = F^ — 1, so this implies 
F < 0. (With F < 0, the argument is similar: then F < 0, so T,j^ > and F = — 1 at the 
center, which implies F < in the neighbourhood.) 

Then suppose that F > and that the limit at ^ of Fi/i? is finite (which means 
that the limit of Fi/Ai^^^ must be finite). As R 00, the contribution from Fi/(xF) 
in (17.21) becomes negligible, and R > Ri eventually. Thus, with F > 0, F must go to 
+00 as -R — > 00, so T,M > 0. In the vicinity of the center, where R < Ri, we must have 
T,M +Fi/{RT) < 0, so Fi < 0. The condition Fi < also follows when F < 0. Thus 
Fi/Ai^^^ < at the center and, by continuity, in some neighbourhood of the center. Then 
we use the same argument as above to show that it leads to a contradiction in every case. 

In the case F < 0, these problems can be avoided. With F < 0, there are two turning 
points, at i? = i?+ and at i? = -R_ > R+. The integral in (17.21) becomes unbounded at 
both and -R_. The signs of the two infinities can then be imposed independently. The 
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argument used after (110. ip still applies, so Fi cannot secure the correct sign at R+. Thus 
we have to assume that the hmit of Fi/ R as M. — > is zero, so that T,_m < secures the 
right behaviour at the inner turning point. Then we require that at -R_ the inequahty 
1 — (c'^/G) {Q,n'^ + QQ,NN ) > — rr^A^ R- holds, which reads exphcitly: 



1 



c 
G 



{Q,n'^ + QQ,NN ) > 



rr, 



M 



2E 
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Since T,m< and E 
necessary condition; 



< in a neighbourhood of = 0, this implies the following 



{Q,N + QQ,NN ) > 



rr. 



M 



2E 



M. 



(10.3) 



Since the integral in (17.21) goes to +oo at both turning points, it must have at least one 
local minimum somewhere in (-R+, i?_). If — cxd < /min < at the smallest of these minima, 
then this can be corrected by the choice of g{M)- Also, / equals to zero in at least one 
point - at -R = -Ri, where R^ < Ri < Thus, / cannot be positive everywhere, so 
necessarily g{Ai) > to avoid F = anywhere. 

In summary, a weakly charged configuration that could bounce singularity-free through 
the R-N wormhole must obey the following necessary conditions in a neighbourhood of 
the center: 

(1) E' < (as E > was eliminated above); 

(2) > — 1/2 (since r is the primary arbitrary function, and E was defined as an 
auxiliary quantity by 2E = r^ — ijj; 

(3) lim^-^rc Fi/ A4^^^ = 0, in consequence of lim^^rc Fi/R = and of (16. 2p (recall that 
/3 7^0 in (O). 

(4) T,_\4 < (to secure the right sign of F at the inner turning point); 

(5) Q,N^ < G/c^ at iV > and Q^n"^ = G/c^ at = (the defining condition for 
weakly charged dust); 

(6) M = — QQ,N r > (a necessary condition for no BB/BC with E' < 0, as shown 
in Sees. 4 and 8); 

(7) — 2EQ'^ {Q-iN^ ~ C!/c*^ > (a necessary condition for the existence of a solution 
of (12.231) . see (14.21) : equality leads to a static solution); 

(8) condition (110. 3p (a necessary condition for (110. 2p to be obeyed); 

(9) condition (110. 2p (to secure the right sign of F at the outer turning point). 
In addition to that, all the regularity conditions of sec. 6 must be obeyed. 

Conditions (1) - (9) must hold in a neighbourhood of the center. At the center, the 
left-hand sides of conditions (1) and (6 - 9) must have zero limits. 

Condition (110. 2p should hold at the state of maximal expansion. However, for those 
shells that contain small mass (are sufficiently close to the center), the value of i?_ will 
be arbitrarily small. Thus (110. 2p must in fact hold down to arbitrarily small values of M., 
with only = being excluded, where both sides of (110. 2p should be zero. 

Now, a glance at eqs. (17.11) - (17.21) suffices to see that the same conditions (1) - (9) 
must apply in the expansion phase, where > and VF < 0, if shell crossings are to 



■^The condition E > —1/2 finds a clearer explanation in the Lcmaitre ~ Tolman (zero charge) limit, 
where it follows from the requirement of the right signature. It is also equivalent to the statement that 
there can be no turning points for radial motion of uncharged dust inside the apparent horizon. 
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be avoided. The two functions g are related to each other by (19.61) . and it may happen 
that if we guaranteed F = + (^c > in collapse, then 1^ + Qe = —h + 9e will refuse to 
be negative in some range during expansion, thus indicating a shell crossing. But if the 
bounce is time-symmetric, so that = —Qc, then, having guaranteed no shell crossings in 
collapse, we know that there will be no shell crossings in the subsequent expansion phase. 
However, the next bounce will in general no longer be time-symmetric, and shell crossings 
will appear in later cycles. The shell crossings might be avoided for ever only if all bounces 
were time-symmetric. Equation (18. 8p implies that such models might exist: the period of 
oscillations T = 27rM/(— 2i?)^/^ will be independent of M when M/(— 2i?)^/^ is constant. 
Together with being constant, such a condition reduces the Lemaitre - Tolman model 
to the Friedmann limit but here we still have the function Q that generates the electric 
field. With T = Tg = const, the function Q{N), in general arbitrary, has to obey 
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(10.4) 



and it can still obey the regularity conditions Q = and \Q,n I = VG/c^ at the center, 
if and F obey their regularity conditions. Also, (110. 4p does not lead to any simple 
contradiction with conditions (5 - 9). An explicit example of a configuration obeying 
(110. 4p (or a proof that it does not exist) remains to be found. 

Equation (110.41) is a consequence of the requirement that the period with respect to 
the proper time s is independent of Ai. It is easy to find from (19.11) that the period with 
respect to the time coordinate t (equal to twice the integral of (19. ip from -Rmin = R+ to 
^max = R-) is 27re-^(^'^(^))/2j\^/(-2E)3/2, where R{M) is an intermediate value of R 
between -Rmin and -Rmax (from the mean value theorem for integrals). Thus, the constancy 
of the period in t also imposes an additional equation on Q. Here, however, an example 
could be found only numerically. 

It should be noted that E need not have the same sign in the whole volume of the 
charged dust. However, the change of sign of E will have no influence on the conditions 
stated above. If it gets positive in a vicinity of the center, then nonsingular bounce cannot 
occur, as shown above. If it gets negative in a vicinity of the center, then it is negative in 
an open range of values of Ai. Then, within the same open set there will be outer turning 
points of some mass-shells, so (110. 3p and (110.21) must anyway hold close to the center. 
However, if E becomes positive at some A4 = Mo, then there are no outer turning points 
at M. > A4q, and the inequahties (I10.2p and (110. Sp need not be imposed in that region. 

We have not yet excluded one more pathological situation that could occur with E < 0. 
The arbitrary functions may be chosen so that the period of oscillations tends to zero as 
Ai ^ 0, as in Fig. [HI The thicker curve (call it Ej^) goes through the minima of the 
evolution curves R{t, A4), each of them corresponding to a fixed value of A4. In the figure, 
the slope a{Ai) of Em is at each point smaller than the slope of the evolution curve that 
passes through the same point, and the shell crossings are inevitable. We could try to 
avoid it by requiring that a{Ai) is sufficiently large, but this would not work: the period 
of a curve with larger Ai would be larger than the period of a curve with smaller Ai, 
and an intersection would be sure to occur after several cycles. With the period tending 
to zero as ^0, the neighbourhood of the center At = would be densely filled with 
shell crossings. Thus, the only way to avoid such a pathology is to choose the arbitrary 
functions so that the period tends to a nonzero value as ^ 0. 
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Figure 6: A possible pathology of an oscillating model with E < 0: the period of oscillations tends to 
zero as the center of symmetry is approached (the curves are graphs of the functions i?(t, A4) for different 
values of A4). The thicker curve goes through the minima of the evolution graphs. (The model shown in 
the figure is time-symmetric with respect to the bounce in the middle of the axis.) 



Equation fl8.8p will help us in choosing the right functions. We assume that nonzero 
intervals of proper time along each worldline are mapped into nonzero intervals of the 
comoving time-coordinate and vice versa (i.e. that the metric function e*^ in does not tend 
to zero or to infinity at the axis = 0). Thus, if we choose the functions so that the 
period in s in eq. (18. 8p is nonzero, then the period in t will also be nonzero. Note then 
that if E obeys (16.41) with a nonzero constant, then, in consequence of condition (5), the 
period of i? in s will tend to zero as ^ 0. To avoid this, E must tend to zero at the 
center faster than A^^/^. 

11 An example 

In order to prove that conditions (1) - (9) from the previous section are not mutually 
contradictory, we shall now provide an example of a model that obeys them. It is simply 
the first model that has successfully passed all the tests, and we do not claim that it is 
physically important or realistic - it is only meant to be a proof of existence. 

Since the functions appearing in the inequalities are rather complicated, the proof that 
the inequalities are all obeyed will be given mainly by numerical graphs. The connection 
between F, N and Al provided by (15.161) causes that even quite simple functions E{N) 
that obey the regularity conditions lead to complicated expressions for Ai{N). However, 
we will verify by exact methods that in a vicinity of the center and at infinity the functions 
indeed behave in the same way as the graphs indicate. 

Choose, in agreement with the regularity conditions and the weak-charge conditions: 
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where q = ±1, to allow for any sign of the charge, x '= N/Nq, Nq is a constant, and 

p{x) = x/il + xf. (11.2) 

Then 



Q'-=^^(^- (11-3) 

This charge is zero at x = 0, its absolute value increases at first, but then decreases and 
tends to zero as x oo. Thus Q,n, which is proportional to the charge density (see 
fl2.20p ). changes sign at x = 1, but it obeys condition (5) everywhere. 
Further, we have 

Fiix) 1 - ^ (Q,^^ + QQ,^^) = 1 - + \ (11.4) 

Fig. [7] shows the shape of the functions p{x), dp/dx and Fi{x) (for all these functions it 
can be easily verified that the graphs show their behaviour faithfully, with no important 
details being hidden beyond the range of the graph or in a small vicinity of zero). 

We choose now the function E{N) so that the period 27rM/(— 2i?)^/^ in (18. 8p has a 
nonzero limit at — 0. In order to keep —1/2 < E < Oin the whole range, and 
to avoid F and Ai = {G/c^) j{l/T)dN being too complicated, we choose the trial form 
2E = —hx""/ (1 + fex'*), with a, h being constant, and we find from (18.81) that the period will 
have a nonzero (actually, infinite) limit at = if a = 5/3. This will also guarantee, via 
(EHD, that the limit of R/ M^'^ at x ^ will be a nonzero function independent of 7W, 
just as (16. 2p requires. Thus 

6x^/^ 

2E = -- —nr. 11.5 

1 + 6x^/3 ^ ' 

With such the limiting period of oscillations in (18. 8p is infinite, and 

F(x) = , \ , (11.6) 

which obeys condition (4), and, further: 

(11,7) 

^ -'0 V 1 + hx'^'^ ^ 

The /i(x) is an increasing function for all x > 0; its graph would almost coincide with the 
graph of F2 shown in Fig. [3 but see the inset (actually, /i(x) > F2{x) for < x < 1, the 
two curves intersect at x = 1, and for x > 1 ^{x) < F2{x)). Note that such obeys 
\\mx^QM./x = GNq/c^ 7^ 0, while lim^^o Fi/x = 12, so condition (3) is obeyed. We have 
now: 

GN 

M = M-QQ,Nr = ^ F2{x), 
(l + x)Vl + &x5/3 
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dp/dx 



2 4 6 8 10 

Figure 7: Graphs of the functions p{x), dp/dx, Fi{x), F2(x) and ^3(2;). Note that -1 < dp/dx < +1 
and F2{x) > 0, ^3(2;) > for all x > 0, thus conditions (5), (6) and (7) are fulfilled. All the functions are 
zero at X = 0. As a; ^ 00, F2 and F3 tend to 00, Fi tends asymptotically to 1, while p and dp/dx tend 
asymptotically to 0. The inset shows the functions p{x), ^2(2^) and fj.{x) in a vicinity of a; = 0. It shows 
that, with our chosen value b — 2.5, fi{x) > p{x) for a; > 0. 



It is easy to verify that F2 > for all x > ojfl so condition (6) is fulfilled. The graph of F2 
is shown in Fig. [7|with b = 2.5 (why this value - see below). 
Condition (7) is equivalent to 



F.(x) > 0, 



where 



F,{x)''=' -2Ep'{l-pJ) 



Fs(x) F2 



^3,11/3 



(x) - F.ix) 



(l + xY 



;ii.9) 



:ii.io) 



The graph of F^lx) with b = 2.5 is shown in Fig. [71 Since we have verified that M > 0, 
E <0 and Q,n'^ < G/c^, equation (111.91) ) is equivalent to F3 > 0, where F^^x) =^ -^2(3^) — 
^yFQ(x). Note that -^3(0) = 0. We divide Fs{x) by and observe that the second term 
in the resulting expression is zero at x = 0, while liniaj^o -^2(a^)/a;^ = 6 > 0. Thus, F3 > 
in a vicinity of x = 0, and, similarly to F2{x), tends to +00 as x — > 00. 



Since T,m = r,Ar N,m, after using (15.161) . we obtain: 



rr. 



M 



2E GNo QxVl + 6x5/3 ■ 



and so condition (8) becomes: 
F4(x) > 0, where 



F,{x) Fi(x) 



6xVl + 6x^/3 



F2{x) 



fll.ll) 



:ii.i2) 



^We have F2{x) > fi{x) -x/ (1 + a;) Vl + fox-Vs 
so F2 > F2 > for all a; > 0. 



dof 



F2 for aU a: > 0, and dF2/da: > for all x > 0, 
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Figure 8: Graphs of the functions F4{x) and Fz{x) with b = 2.5. The graphs of Fi{x), F2{x) and ^^3(2;) 
are shown for comparison and for scale. Both F^i^x) and ^5(2;) asymptoticaUy tend to 1 as x — > cx) and 
both have vertical tangents at a; = 0. 



=0 



We have -^4(0) = 0. Knowing that hm^^^o F2/X = 6, we easily calculate that dF4/dx|^^ 
7 > 0. Thus F^lx) is increasing in a vicinity of x = 0, i.e. F^lx) > at least for some 
X > 0. Also, Fi{x) — > 1. The graph of F^lx) with b = 2.5 is shown in Fig. [HI 

a;— >oo 

Finally, condition (9), i.e. eq. (110.21) . can be written as 



dof 



F,{x) > 0, where F,{x) = F^{x) ^==j= H^) + VFsix) , (11.13) 

oxv 1 + bx^''^ L -I 

(see eq. (111.101) for the definition of Fq{x)). Since we have already verified that (I11.12p 
holds, we can take the term with the square root in (111.131) to the r.h.s. in F5 > and 
then square both sides of the inequality. The result can be rewritten as: 

/ 5 \ 25 

We note that lima;^o -^i/a; = 12, lima;^o -^2/2;^ = 6 and limx^o F^/x^^^'^ = 8b. Then we 
divide (111.141) by x"^ and take the limit of the resulting expression at x ^ 0. The result 
is 24 > 0. Hence, (111.141) is obeyed in a vicinity of x = 0, and, consequently, F5 > is 
obeyed in that vicinity, too. It is easy to verify that F^lx) — > 1. The graph of -F5 with 

a;— »oo 

6 = 2.5 is shown in Fig. [HI 

The reason for choosing b = 2.5 was this: the graph of F5 is sensitive to the value of b. 
With decreasing b, the local minimum of F5 becomes smaller, and for b small enough (for 
example, 6 = 0.75) F5 < around the minimum. With b > 2.5, the minimum is clearly 
positive. 

Graphs do not constitute a definitive proof that a function is positive in the whole 
infinite range. Fine details can be hidden beyond the range of the figures or at a very 
small scale around the points where the functions approach zero. In those regions the 
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0.5 1 1.5 2 2.5 3 3.5 4 

Figure 9; The maximal radius achieved in each cycle, and the radius of the outer R-N horizon, 
Rh+, as functions of mass. The numbers on the axes are multiples of GNo/c^. At this scale, the radius of 
the inner horizon and the minimal radius seem to coincide with the a;-axis, but see Fig. 1101 As required 
by the equations of motion, we have i?_ > Rh+ > Rh- > R+ for all masses. Note that i?_ and Rh+ are 
increasing all the way. 



functions can go below zero, while we cannot see it. However, we have verified by exact 
methods that there exists a neighbourhood ?7 of x = in which the functions Fi, . . . , F5 
are positive for x > 0. Thus, even if some of the functions become negative outside U, we 
can cut away a finite ball of the charged dust with the radius xo smaller than the radius 
of U and match the charged dust to the Reissner - Nordstrom solution at x = xq. In this 
way, we will obtain a finite charged body of dust that can go through a minimal size and 
bounce without encountering any singularity. 

It can be easily verified that in our example {G/c^)Q'^ < M.^ for all x > provided 
h < 25.3 (hint: p,x< is equivalent to [(1 +x)^/(l - x)^ - 1]^^^ > 53/2^5/2 _ ^^le 
inset in Fig. [7]). Thus, if a finite sphere is cut out of this configuration and matched to the 
Reissner - Nordstrom solution, the exterior R-N metric will have < m^, and horizons 
will exist in itj§ As is well-known from the studies of the Reissner-Nordstrom metric [S], 
the surface of a charged dust sphere matched to the R-N metric moves according to the 
same equation as a charged particle moving radially in an R-N spacetime. Consequently, 
the reversal of collapse to expansion cannot occur between the two R-N horizons. Our 
example is consistent with this. Fig. M shows the radii of the outer horizon and of the 
maximal radius achieved by a given Ai shell as functions of x = N/Nq. At the scale of 
Fig. [9l the minimal radius and the radius of the inner horizon seem to coincide with each 
other and with the x-axis. Fig. [10] shows a closeup view of the two curves, and shows 
that indeed the inner horizon has a smaller radius everywhere, except at the maximum, 
where the two curves are tangent to each other. This is no numerical coincidence - the 
two curves must be tangent at every local extremum of R+. 

The limit of the period (in s) of oscillations in (18.81) at x — > is infinite. Thus, the 
limiting period in the time coordinate will also be infinite. Fig. [11] shows the numerically 

^But GQ^ /{c'^Ai^) —> 1 as X ^ 0, so as the sphere becomes smaller, the outer R-N metric tends to the 
critical one, = w? . 
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Figure 10: The radius of the inner horizon, Rh- (upper curve), and the minimal radius, R+ (lower 
curve), as functions of mass. The numbers on the axes are multiples of GNq/c^. Both radii first increase 
with mass, but then decrease. The figure shows that i?+ < Rh- everywhere except at the maximum, 
where the two curves are tangent to each other. The inset is a closeup view of the neighbourhood of the 
maximum. 



calculated period in t as a function of A^, in the range x G [0, 4) (for x > 4 the calculation 
becomes progressively more difficult). As the figure shows, the period decreases with 
increasing x in this whole range, so we will assume x = 4.0 as a "practical" point of 
minimum of T. A real minimum must occur somewhere at x > 4.0 because the period 
increases to +oo at both x — > and x ^ +oo. 

Finally, Fig. [T^ shows a collection of curves R{Ai,t) corresponding to different values 
of A4, calculated numerically by solving the set fl2.2ip - (12.231) . The configuration is time- 
symmetric with respect to the instant t = 0. In agreement with Figs. [9|, [10] and [11], the 
maximal and minimal radii achieved in each cycle are increasing functions of mass, while 
the period is a decreasing function^ 

Figs. [12] and [TT] are drawn on the basis of numerical solutions of the set (12.221) - (12.231) . 
As (12.221) shows, for small masses, where Q ~ 0, C is nearly constant, i.e. nearly zero 
(because, as mentioned earlier, t can be chosen so that C = at the center). Thus, the 
dependence of i? on t for small values of x should be very similar to the dependence of 
R on s. Comparison of the corresponding graphs confirms this: the graphs of R{A4,s) 
look identical to those of Fig. [T2] and the period as a function of s has a graph that looks 
identical to Fig. [TT] 

As predicted in Sec. [TDl each mass shell avoids shell crossings throughout the first 
expansion phase after the time-symmetric bounce, but then experiences crossings after 

^The range of masses in the figure is x G [0.01,0.1], where Figs. [91 [TUl and [Til show that the radii and 
the period should really behave the way they do in Fig. [T^l For Figs. [TT] and [T^ the assumed value of b is 
0.001. 
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Figure 11: The period (in the time coordinate) as a function of a; = N/Nq. 



4500 



tr. 




20 40 60 80 100 120 140 

t 



Figure 12: The curves R{A4,t) corresponding to several vahies of A4. The mass increases uniformly 
from X = 0.01 on the lowest curve to x — 0.1 on the highest curve. The range of masses was chosen so 
as to make the figure readable; the real range of periods and of amplitudes can be very large. See more 
explanation in the text. The bounce is always smooth and at a nonzero value of R. (The figure suggests 
otherwise, but this is only an illusion created by the scale.) 
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going through the maximal size. In order to avoid shell crossings in the whole volume for 
the whole expansion phase of the outermost shell, the radius of the dust ball cannot be 
too large. If it is very large, then the period of oscillations of the outermost shells will 
also become very large. Then, the time by which first shell crossings appear inside the ball 
will become a small fraction of the duration of the expansion phase of the outer surface, 
i.e. shell crossings will appear before the surface of the ball emerges from inside the outer 
horizon. 

The evolution of our configuration is summarised in the Penrose diagram in Fig. [T3l 
The diagram is written into the background of the Penrose diagram for the maximally 
extended Reissner-Nordstrom spacetime (thin lines). C is the center of symmetry, Sb is 
the surface of the charged ball, Srn is the Reissner-Nordstrom singularity. The interior 
of the body is encompassed by the lines C, E, Sb and B; no singularity occurs within this 
area. Lines B and E connect the points in spacetime where the shell crossings occur at 
different mass shells. Nl (N2) are the past- (future-) directed null geodesies emanating 
from the points in which the shell crossings reach the surface of the body (compare Fig. 
[T2|) . The line Sb should be identified with the uppermost curve in Fig. [121 The top end 
of Sb is where the corresponding curve in Fig. [12] first crosses another curve, the middle 
point of Sb is at t = in Fig. [T2j 

It would be interesting to have an example of a configuration that can pulsate for ever, 
avoiding shell crossings in all of its collapse/expansion phases. However, our example does 
not obey (110.41) - substitution of our functions into (110. 4p leads to a clearcut contradiction. 
Thus, it is impossible to impose on our example the condition that the period (in the 
proper time s) of pusations is independent of Ai. Whether such a permanently pulsating 
singularity-free configuration exists at all is a problem to be investigated in the future. 

We recall that our example was only meant to demonstrate that the inequalities in 
conditions (1) - (9) are not mutually contradictory and allow a solution. Since it turned 
out that the set of models is not empty, more examples should exist, and now it is a 
challenge to explore other possibihties. 

In the Friedmann limit (i.e. zero charge density, zero central charge and homogeneous 
mass density), a model with E < goes over to the k > Friedmann model, for which the 
most natural topology of spatial sections is that of S^. Here we have a model that, if it 
were not matched to an R-N spacetime, would have spatial sections of infinite volume (as 
attested by the fact that it can contain an infinite amount of rest mass). Lemaitre-Tolman 
models with such spatial topology are known and understood geometrically [15] . However, 
the limiting transition from such a model to the A; > Friedmann model can be done only 
locally, i.e. in finite volumes, and involves a discontinuity in the arbitrary functions. Thus, 
it would be desirable to find a charged dust model with E < that would have complete 
spatial sections of finite volume. Curiously, it turned out to be much more difficult in this 
case to fulfil conditions (7) - (9), and no example of such a model has been found. (The 
authors do not wish to imply that such an example does not exist, this is simply a problem 
to be solved.) 
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Figure 13: A schematic Penrose diagram for the configuration defined by eqs. (jll.ip - (|11.2|) and ()11.5|) . 
See explanation in the text. 
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12 Can such an object exist in the real Universe? 



There is a belief in the astronomical community that ordinary astronomical objects, such as 
stars and galaxies, have zero net charge. Quantitative estimates based on measurements 
and observations are, however, hard to come by. The best that can be found in the 
literature are theoretical considerations on how charges could be separated within the body 
of a star and what the maximal charge could be. Before we compare our results with those 
predictions, we shall first choose the most favourable parameters for our configuration. 

It is clear that the smaller the net charge, the greater the chance that such an object 
might exist. For our object, the absolute value of the charge first increases from zero in the 
center to the maximum \/GNq/ {Ac^) achieved at x = 1, and then keeps decreasing all the 
way to zero as x ^ oo. It needs to be explained how such a charge could be concentrated 
toward the center. We shall not dwell on this mechanism because stars are obviously 
not composed of dust. In our case the concentration is a dynamical effect, achieved by 
fine-tuning of the initial positions and velocities. 

As seen from Fig. [11], the period of oscillations decreases with mass, reaching a mini- 
mum at X > 4.0. Then, Fig. [9] shows that at x = 4.0, the maximal radius is slightly larger 
than the radius of the outer horizon. The period in s, calculated from (18. Sp . 1t\:M j i^—IEf'^'^ ^ 
goes to oo both at x — and at x — > cxd, going through a minimum at x ~ 1.0. The 
period in t should behave similarly, i.e. will start increasing at some x = Xm- We cannot 
take the radius of our object larger than Xm because if it is large, then the period will 
be large, and shell crossings will appear inside the object before its surface emerges from 
the outer horizon. This would be unsatisfactory - we want the object to distinctly emerge 
from the horizon before any singularities destroy it. Thus, the largest total mass that we 
can assume corresponds to x ~ 4.0. At x = 4.0, the charge is \Q\ = 0.16\/GNo/c^, while 
the active mass, in physical unitsjfl is c^AA/G = 1.68No/c^. The ratio of charge to mass 
is thus G\Q\/{c'^A4) = 0.095-\/G in electrostatic units (in which the unit of charge is y/g 
cm^/^/s). This makes 0.82 x 10~^^ coulombs per gram. For the whole Earth, this ratio is 
0.502 X 10~^^ C/g [16j. Compared to the charge of the Earth, the charge of our object is 
thus enormous. However, for a neutron star of 1 solar mass, the authors of Ref. pT] found 
that the total charge might be 10^° C, which makes 5.03 x 10~^^ C/g - distinctly more 
than in our object. Thus, the possibility to find a real object with charge and mass similar 
to our example is not outlandish]^ Only the mechanism of charge separation remains a 
problem. 

We stress again that the functions defining our object were chosen rather accidentally, 
and lower ratios of charge to mass at the surface can possibly be achieved. But the limiting 
ratio of charge to mass at the center must be ^/G, which makes 8.616 x 10^^^ C/g. 

The period of oscillation, found from (18. 8p . is T = 27rM/ (— 2i?)^/^, but this is expressed 
in geometrical units. In physical units, the period is T/c. Assuming that the surface of 
the object is at x = 4.0, taking the corresponding mass to be = 1.68GA^o/c^ taking 
b = 2.5, we calculate 

T/c = 2.78 X 10"^^ (No/c^) , 
which, for solar mass Nq/c^ = 1.989 x lO^^g, gives 5.52 x lO^^s - the same order of 
^What we have been calling 'active mass' M is in reality {G/c^jx (mass). 

^Note that we have not fixed the mass of our object - the constant A^o is still arbitrary and can be 
made whatever we wish. 



33 



magnitude as the time of collapse to the Schwarzschild singularity by a neutral star of one 
solar mass that has crossed the horizon. Needless to recall, for an observer at infinity in 
the R-N spacetime, the time it takes a collapsing object to reach its horizon is infinite. 

13 Conclusions and possible further research 

We have shown that it is possible to set up such initial conditions for a charged dust 
sphere of a finite radius that its outer surface completes one full cycle of pulsation, while 
no singularity appears either at the surface or anywhere inside it. The sequence of events 
is this: 

1. The initial instant is when the outer surface of the sphere is close to the maximum 
expansion. (When it is too late after the maximum, shell crossings in the interior will 
appear before the full cycle is completed; when it is too early before the maximum, shell 
crossings in the interior will exist already at the initial instant). Other surfaces of constant 
mass reach their maxima of expansion at (slightly) different instants of the time coordinate 
t. The outer surface at the initial instant is outside the outer Reissner-Nordstrom horizon. 

2. During collapse, the outer surface plunges first through the outer R-N horizon, 
then through the inner R-N horizon, and bounces at a nonzero radius. No shell crossings 
appear anywhere inside the sphere during the collapse phase, and none of the constant-mass 
spheres collapse to zero radius, i.e. there are no singularities in this phase. 

3. For all constant-mass shells, the bounce at minimal size is simultaneous in the time- 
coordinate t. This means that the evolution is time-symmetric with respect to this instant. 
This makes the problem technically simpler: if there were no singularities during collapse, 
there will be none during the next expansion phase. 

4. The expanding outer surface goes out of the inner R-N horizon, then out of the 
outer R-N horizon and expands up to the maximal radius. As follows from the analysis of 
the maximally extended R-N spacetime, after emerging from the outer R-N horizon, the 
surface is in a different asymptotically flat sheet of the extended R-N spacetime. 

5. After the outer surface goes through the maximal radius, during the second collapse 
phase, shell crossing singularities appear inside the ball and its further evolution cannot 
be followed. 

With the most favourable value of mass, the ratio of total charge to total mass of our 
dust ball is only slightly larger than the corresponding theoretically estimated maximum 
for charged neutron stars (see our Sec. 12 and Ref. [H]). 

This is the description of our chosen example. However, other choices of the arbitrary 
functions are possible that might improve some of the characteristics of our model. For 
example, it would be desirable to avoid the matching to the R-N spacetime, so that the 
dust distribution can extend over the whole space (infinite or closed) - then the solution 
could be investigated as a possible model of the Universe with a localised charged object 
in it. Another desirable generalisation would be to make all bounces time- symmetric, 
so that the model oscillates singularity-free for ever. We shall now briefly discuss these 
possibilities. Recall that one of the necessary conditions for the absence of singularities is 
E < 0, and we shall assume this in the following. 

Let us recall the necessary conditions for avoiding the BB/BC and shell-crossing sin- 
gularities. Two conditions are most important: 
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1. The absolute value of the ratio of charge density to mass density at the center 
must be (in the units used here) \/G/c (see (12.201) ). Otherwise, either shell crossings are 
inevitable or the ratio of charge density to mass density must be large throughout the 
volume to prevent BC. 

2. The limit of the period of oscillations at zero radius must not be zero, so that the 
pathological situation of Fig. [6] does not occur. 

In addition to that, the following conditions must be obeyed: 

3. The regularity conditions at the center of Sec. O which are the required values of 
various functions at = 0. 

4. The 9 conditions listed in Sec. 10, most of which are inequalities to hold throughout 
a neighbourhood of the center. Of these, the functions in conditions (1) and (6 - 9) must 
have zero limits at the center. 

All these conditions are difficult to fulfil because there are only two arbitrary functions 
in the model: (z) the charge distribution Q{N), and (ii) the energy distribution E{N) 
(which can equivalently be defined by specifying r(A^) or A^(A^)). Our choice here was 
the first {Q, E) pair, found by trial and error, that obeys all the requirements, but only in 
a limited volume and in a limited time interval. 

To avoid the artificially limited volume, one might try to match the model to a space- 
time different from R-N. For example, it could be the Vickers model with the charge 
density becoming strictly zero at a certain distance from the center. Then, in the outer 
region, Q,n= 0, which is sufficient to prevent the BB/BC singularity. Since the main 
difficulty in avoiding shell crossings was close to the center, and the Q,n = region would 
not extend to the center, chances are that shell crossings could be avoided as well. In the 
extreme case, the outer region could be the Lemaitre-Tolman spacetime. Since the L-T 
spacetime has strictly zero electric field, the function Q{N) would have to be zero at the 
surface of the charged region. The zero total charge would be an advantage from the point 
of view of astronomy. We do not know if such a configuration exists. 

To avoid the limited time interval, the bounce at minimal radius should be simultaneous 
in t for all constant-mass shells. This means that the period of oscillations (measured by 
the time-coordinate t) would have to be independent of mass. Technically, this means that 
the integral of (19. ip from the minimal i?+ to the maximal i?_ (both given by (14. 3p ) should 
be independent of Ai. We have not been able to tackle this condition. 

The condition to avoid zero limit of the oscillation period at the center is, however, 
quite simple. It is reasonable to require that the function C has a finite value at the center, 
where = - without this assumption, finite intervals of the time coordinate t would 
correspond to infinite intervals of the proper time s. If this is the case, then it is enough 
to choose such functions that the period in s has a nonzero limit at A^ = 0. From (18.80 . 
the condition for this is that M / {—2E)^^'^ has a nonzero limit at A^ ^ 0. As noted at the 
very end of Sec. 10, the necessary condition for this is: the limit calculated in (16. 4p is zero. 
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